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ABSTRACT 



Cool stars are known to produce flares probably as a result of magnetic reconnection in their outer atmospheres. We present simulta- 
neous XMM -Newton optical V band and X-ray observations of the M8 dwarf LP 4 1 2-3 1 . During the observation a giant flare occurred, 
with an optical amplitude of 6 mag and total energy of 3 ■ 10 32 erg in each the V band and soft X-rays. Both flare onset and flare decay 
were completely covered in both wavebands with a temporal resolution of 20 s, allowing for a determination of the flare time scales as 
well as a study of the temperature evolution of the flaring plasma. The data are consistent with an impulsive energy release followed 
by radiative cooling without further energy release during the decay phase. Our analysis suggests that the optical flare originates from 
a small fraction of the surface of LP412-31, while the characteristic scale size of the flaring X-ray plasma is of the order of the stellar 
radius or larger. The absence of any small-scale variability in the light curve suggests a non-standard flare number energy distribution. 

Key words. X-rays: stars, stars: activity, flare, low-mass, brown dwarfs, coronae, stars: individual: LP412-31 



1. Introduction 

Flares, i.e., short-term enhancements of the emission in various 
wavebands, belong to the most obvious manifestations of solar 
and stellar magnetic activity. During a solar white-light flare the 
footpoints of magnetic loops light up in optical and ultravio- 
let emission lines and continuum. Intense X-ray emission goes 
along with such events, and is ascribed to the filling of these 
loops with heated plasma evaporated off the chromospheric lay- 
ers along the magnetic field lines into the corona (the 'chro - 
mospheric evaporation scenario', e.g.. lAntonucci et all I1984). 
Depending on the flare frequency distribution quiescent high- 
energy emission might be due to a super-position of unresolv ed 
micro- flares (iKrucker & BenJl998l lAschwan den"e"t all2 000). 

On dM(e) stars ("flare stars") flares are far more frequently 
observed (mostly in U band or Ha) than on the Sun. However, 
because of significant logistic and observational difficulties, si- 
multaneous data taken in different wavebands, indispensable 
to test the validity of the above scenario, has rarely been se- 
cured. The situation becomes even worse for the low-mass end 
of the main-sequence. The fraction of Ha emitters among the 
M dwarf population decreases quickly for ultracool dwarfs be- 
yond spectral type m M7, yet objects of late-M or L spec- 
tral type can produce spectacular Ha flares (Liebert et alJll999i 
Lieb ert et aflboOfl iFuhrmeister & Schmittll2004l) . X-rav emis- 
sion has been observed only from a rather small number of ul- 
tracool dwarfs, and in most of these cases X-ray detection s were 
made only during flares (Fleming et al. 2000; Rutledge et al. 
2000; ISchmitt & Liefkd l2002t IStelzerl l2004l) . None of these 
observations was accompanied by concurrent observations in 
other wavebands. Optical broad-band observations of (ultra- 
cool) dwarfs are more abundant. However, most of the avail- 
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able data do not have the temporal resolution required to fol- 
low the rapid evolution of the photospheric and chromospheric 
emission. The observed magnitude increases in the optical wave- 
bands can be enormous. Some - rare - examples of huge opti- 
cal events on M dwarfs with amplitudes > 6 mag in optical and 
UV bands, respectively, are presented by Pag ano et alJ Jl997h 
an d lRockenfeller et alJEOOoT) . Observational records of nine su- 
perflar es in various wavebands were collected by Schaefe ret alJ 
(2000), but these observations lack any information on time evo- 
lution. 

In this letter we present a giant X-ray flare on the M8 
dwarf LP 4 12-31, recorded with the XMM-Newton satellite with 
a time resolution of a few seconds simultaneously in soft X- 
rays and in the V band. LP 4 12-31 was o riginally identified 
in a proper motion survey of Luvten ( 1979) an d first classi- 
fied as a very late M dwarf by Kirknatric k et alJ (119951) . Since 
then it has frequently been used as reference star for spec- 
tral type M8, although little is actually known about its stel- 
lar par ameters. LP 4 12-31 is located at a di s tance of 14.6 + 
0.1 pc ( iReid & Cruzll2002l) . iMohantv & Basril ( 120031) derive an 
effective temperature of 2550 K, and from the J band mag- 
nitude we obtain Lboi = 2.09 • 10 30 ergs~ 1 usin g the bolo- 
metric correction (B.C.)j = 1.95 (Re id et al]l2001l). A ccording 
to the evolutionary models of Chabrier et al. (2000) LP 41 2- 
3 1 is close to the bottom of the main-sequence (0.04 -0.1 M ) 
and moderately young (« 0.3 Gyr). LP 412-31 belongs to the 
strongest Ha emitters of its spectral type. Its Ha equivalent 
width has been measured at various epochs with values between 
20...80 A. Martin et alJl996llGizis et alJ2000tlReid et alJ2 002: 
Mohantv & B asril2003l) . The magnetic field strength of LP 412- 
31 (> 3.9 kG) is among the highe s t in a sample of ultracool 
dwarfs studied by Reiners & Basri (2006). No reports on X- 
ray emission from LP 412-31 exist in the literature, and in par- 
ticular the ROSAT All-Sky Survey data implies a flux limit of 
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< 2 x 10 13 erg cm 2 s 1 , corresponding to a soft X-ray luminos- 
ity of < 5 ■ 10 27 ergs~'. 

2. Observations and data analysis 

We observed LP 4 12-31 for 40ksec with XMM-Newton 1 
(obs-id0300170L01), and analysed the data with the XMM- 
Newton Science Analysis system (SAS) version 6.5.0. 

The thin filter was inserted for all three European Photon 
Imaging Camera (EPIC) instruments. In this paper we use only 
data from the EPIC/pn CCD. Data processing included standard 
filtering of the event list. After removing time intervals of high 
background using a screening algorithm, the EPIC/pn yielded 
25ksec of useful science data. Visual inspection of the X-ray 
image showed an unexpectedly bright X-ray source coincident 
with LP 4 12-31. For the temporal and spectral analysis of this 
source we used the source photons extracted from a 30" cir- 
cle centered on the X-ray position of LP 412-3 1 and background 
photons from a circle of equal size positioned on a source-free 
adjacent position. 

The Optical Monitor (OM) onboard XMM-Newton was oper- 
ated in the standard Image Fast Mode with the V band filter. For 
this work we use only data from the central fast-mode window 
of 10.5" x 10.5" pixels read out with a time resolution 20 s. For 
technical reasons the OM observation had to be divided into nine 
exposures (7 x 4400 s, 1 x 3700 s, and 1 x 2960 s) separated by 
an overhead of ~ 0.3ksec. Thus, there is almost complete time 
coverage of the X-ray observation in the V band with high time 
resolution. The OM data from the central fast-mode window was 
processed with standard SAS procedures, resulting in a sequence 
of V band count rates with time resolution of 20 s. For the con- 
version of OM count rates into magnitudes and fluxes we used 
the numbers given in the XMM-Newton SAS Users' Guide 2 . 

An inspection of the lightcurves of both the X-ray and opti- 
cal data showed the presence of a huge flare at about one third 
into the observation, but a weak signal without any signs of vari- 
ability during the rest of the observation. Substantial precursor 
activity before the actual flare onset is clearly visible in the X-ray 
data (EPIC/pn count rate enhanced by a factor m 30 with respect 
to the pre- and post-flare quiescent emission). Fig. ^ shows a 
fraction of the recorded lightcurves, starting with the X-ray flare 
precursor, during which no significant activity is seen in the OM 
data. In Table^we summarize count rates, magnitudes, and lu- 
minosities measured during the different observed activity states. 

For the quiescent state before and after the flare (not shown 
in Fig. [Q, we measure an X-ray luminosity of logL X; q ~ 
27.2 erg s , corresponding to a fractional X-ray luminosity of 
log (Lx,q/^boi) 3.1. This represents the first detection of (qui- 
escent) X-ray emission for LP 4 12-31. Our measurement puts 
LP 412-31 at the upper end of the activity range for flare stars, 
not unexpectedly considering its strong and variable Ha emis- 
sion and strong magnetic field. The quiescent OM V band count 
rate translates into a V band magn itude of ~ 19 ma g, in good 
agreement with the literature value (Dahn et al. 2002). 

Both visual inspection of the lightcurves and Spearman's 
and Kendall's rank correlation tests are consistent with a time 
delay of < 20 - 40 s (1 - 2 bins) between the optical and X- 
ray bands. Both EPIC/pn and OM count rates increase by fac- 
tors of 200 - 300 from quiescence to their peak luminosities of 

1 The satellite and its instruments have been described in a Special 
Issue of Astronomy & Astrophysics Vol. 365 (Jan. 2001) 

2 http://xmm.vilspa.esa.es/external/xmm_user_support/documentation 
/uhb/index.html 
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Fig. 1. Fraction of the XMM-Newton lightcurves of LP 412- 
31 showing the giant flare and the precursor state. Both EPIC/pn 
and OM data are binned to the 20 s time resolution given by the 
OM. The EPIC/pn lightcurve is overlaid by the blocks derived 
in the ML analysis. The last block is not seen because of the 
restricted x-scale. 

4.7 ■ 10 29 erg s _1 and 2.8 • 10 30 erg s _1 , respectively. Specifically, 
the V magnitude rises within 40 s to a maximum of 13.1 mag 
(AV sa 6 mag), whence it decays exponentially back to quies- 
cence with decay time scale tom ~ 90 s. The decay of the X-ray 
light curve is far more prolonged. It is fast initially, but slows 
down soon after and can be roughly described by two exponen- 
tials, a first one with a decay time scale similar to the optical 
decay (Ti pn ~ 150 s), and a second exponential that decays with 
T2,p„~ 1000s. 

We analysed the EPIC/pn light curve with a binning inde- 
pe ndent maximum like lihood (ML) method described in detail 
bv lStelzeret all E006). We use this technique first to find pos- 
sible background variations and then subtract the background 
from the events listed at the source position, making use of the 
previously determined segmentation of the background signal. 
The result is a background subtracted net source photon events 
list, to which the ML algorithm is applied. In Fig.^the binned 
lightcurve is overlaid by the constant rate blocks resulting from 
the ML analysis with an imposed minimum number of 500 cts 
per block and a confidence level for change points of 99.9%. 
Since the data before the flare precursor (not shown in Fig. [0 
consists of fewer than 500 cts, we start the ML analysis with 
the flare precursor. This procedure results in a division of the 
lightcurve into ten intervals of constant signal which we use for 
spectral analysis. 

For each of the ten blocks as well as for the first 6 ksec of 
the observation (representing the quiescent emission before the 
flare) a source and background spectrum, a response matrix and 
an ancilliary response file were extracted with standard XMM- 
Newton SAS tools. We first fitted the quiescent spectrum within 
the XSPEC vll.3 environment using a two-component APEC 
thermal spectrum with a photo-absorption term. The spectral fit- 
ting shows the absorbing column to be very low as expected for 
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Table 1. Count rates, magnitudes and luminosities for LP412- 
31 in different activity states. 







pre-flare 1 "' flare peak 


post-flare* 6 ' 


OM 


C 


[cts -1 ] 


0.50 ±0.20 84.83 ±3.47 


0.30 ± 0.08 


V 


[mag] 


18.7 ±0.4 13.1 ±0.1 


19.3 ±0.2 


L V 


[ergs -1 ] 


(1.6 ± 0.5) ■ 10 2S (2.8 ±0.1)- 10 30 


(9.8 ± 2.0) ■ 10 27 


EPIC/pn 


C 


[cts -1 ] 


1.089 ±0.067 8. 100 ±0.636 


0.042 ± 0.003 


L X 


[ergs -1 ] 


(6.1 ± 0.4) ■ 10 28 (4.6 ± 0.4) ■ 10 29 


(1.6 ±0.2)- 10 27 



(<I) defined by ML blocks No. 1 and 2 representing the flare precursor 
activity (enhanced count rate in EPIC/pn with respect to quiescence; see text). 
™ defined by ML block No. 10 comprising all data after the flare decay. 
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Fig. 2. EPIC/pn spectrum and best fit model for flare precursor 
(Blk2; thick lines), flare peak (Blk3) and the end of the flare 
decay(Blk8). 

a nearby star. For the subsequent analysis of the flare spectra we 
froze all parameters of the quiescent model and added a spectral 
model describing the additional flare emission. 

For the flare component, we tested a series of one- 
temperature (1-T) and two-temperature (2-T) thermal models 
with global abundance Z set to different values (Z = 0.3 Z Q , 
Z = 1.0 Z , and Z a free fit variable). The 1-T models turn out to 
be inacceptable for most time segments. The only approach that 
yields a statistically acceptable description and well-constrained 
parameters for all spectra is the 2-T model with Z = 1.0 Z , con- 
sidered the best fit model (see Fig.[5]for a selection of spectra). 
A mean temperature T was computed for each spectrum from 
an emission measure {EM) weighted sum of the two fit compo- 
nents. The total EM is given by the sum of the EMs of both com- 
ponents. The resulting T vs. VEM diagram is plotted in Fig. [21 
and will be used as a diagnostic for the flare evolution. 



3. Discussion and Conclusions 

3. 1 . Optical and X-ray properties of the giant flare 

To derive quantitative information from the optical signal we 
choose a description of the flare photosphere in terms of black- 
body emission with effective temperature as a free parameter. 
Although the flaring photosphere is likely not a perfect black- 
body, at least in the V band there are no strong emission lines 
Jde Jager et al.lfl989ft . and blackbody models have provided a 
reasonable approximatio n for flare continua on earlier M dwarfs 
(Hawlev & Fisher 1992). Using the V filter response curve, the 
luminosity derived from the OM count rate (see Tabled can be 
corrected for the temperature dependent out-of-band flux ("bolo- 
metric correction"). Assuming a photospheric 'flare tempera- 
ture' r e ff = 9500 K, we obtain an optical flare peak luminos- 
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Fig. 3. Time evolution of coronal temperature and emission mea- 
sure. Open circles denote the flare precursor, filled circles the 
flare decay. The dotted line represents a linear least squares fit 
to the flare decay, that yields a slope £ = 1.25. Numbers in the 
plot refer to the constant rate blocks determined with the ML 
analysis. 

ity of L bo i, peak = 4.7 ■ 10 30 ergs -1 ; raising T eff to 15000K in- 
creases this luminosity by a factor of two. With an assumed stel- 
lar radius of 0.12 R Q , the area of the flaring region on the star is 
(3.4 - 13) 10 18 cm 2 , i.e., the flare occupied < 1 % of the surface 
of LP 412-31. This com pares well with the flare area sizes deter- 
mined forJJVCetjde Jager et al. 1989) and for solar white light 
flares dXu et all2006l) . 

We determined the physical parameters of the X-ray emit- 
ting region making use of the hydrodynamic flare model and 
its calibration to the EPIC/pn spectral response described by 
Real e et all (|2p04). This approach implies that we interpret the 
flare on LP 4 12-31 as a compact loop flare. From the measured 
slope of £ = 1 .25 of the flare evolution in the log T - log y/EM 
plane (see Fig. [3} and the observed log r pea k [K] = 7.65, we find 
a (half) loop length of L « 1.4 ■ 10 10 cm, i.e., a value larger than 
the radius of LP 412-31. Assuming that the cross section of the 
flaring loop be given by the size of the optical flaring region, we 
compute a volume of Vn^-e ~ (0.9 - 3.6) • 10 29 cm 3 and - from 
the observed peak EM of 3.4 ■ 10 52 cirT 3 - a mean plasma den- 
sity of (3 - 6) • 10 11 cm -3 . In principle, the expectation of high 
density could be checked with the XMM-Newton RGS data, but 
unfortunately due to the high flare plasma temperature the sig- 
nal in the O VII triplet (which would be sensitive to densities 
in that range) is so weak that no quantitative conclusions can be 
drawn from that data. At any rate, however, the size of the flaring 
region in X-rays is quite large and of the order of the stellar di- 
mensions. Therefore considerable heating also of the non-flaring 
atmosphere must occur as a result of its illumination by the hot 
flaring plasma. Further, the derived length scale L must be filled 
with plasma quite rapidly. Assuming ad hoc expansion veloci- 
ties of 1000 km/s, the typical filling time scales are of the order 
of 50-100 s, which is comparable to the possible observed delay 
between flare onset and flare maximum. Also, considerable line 
shifts and broadening are expected which, however, cannot be 
diagnosed because of low signal-to-noise of the presently avail- 
able data. Finally, we can estimate the total thermal energy con- 
tent £ t h,x = 3nkTV of the X-ray plasma at flare peak. Using the 
above numbers we find /ith.x = (1 — 2) ■ 10 33 erg. 
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In terms of optical amplitude this flare on LP 4 12-31 is 
among the largest recorded in detail throughout the history of 
stellar flare observations and it belongs to the very few flares 
allowing an assessment of the overall energy budget both in 
X-ray and optical light. The emitted flare energy amounts to 
3 • 10 32 erg, each in the V band and in soft X-rays (see Tabled, 
values clearly at the upper end of the range ob served for field 
flare stars but lower than those for T Tauri stars ( Fernande z et alJ 
2004; IStelzer et alJl2006T) . The total radiated energy is thus at 
least E m i x 6 ■ 10 32 erg, and this is obviously a lower limit, since 
energy may be radiated in other wavebands. We further find that 
the energy radiated in the soft X-ray band is roughly of the same 
order as the value derived for the thermal energy content of the 
coronal plasma at flare peak, and hence there appears to be no 
substantial energy release required in the flare decay phase. This 
is consistent with the fact that the observed value for the slope in 
the T vs. yEM diagram is close to the limiting value of £ = 1.3, 
above which t he heating tim e scale is zero in the hydrodynamic 
model of lReale et all fe004l> . 

In the standard flare scenario it is assumed that the total 
energy released in the flare is ultimately derived from energy 
stored in the magnetic field and converted into other forms of 
energy during the flare. Comparing the magnetic energy density 
to the radiated flare energy and flare volume derived above, we 
can compute the mean magnetic field that must be annihilated 
in the corona in order to provide the observed energy. We find 
Banniii ~ (140 - 280) G. The minimal field necessary for confine- 
ment is the equipartition field, B eq as 3 10- 430 G, and obviously 
this field must not be destroyed in the course of the flare progess. 
The total coronal magnetic field prior to flare onset must thus be 
in the range 450-7 10 Gin a volume of V flare * (9-36)- 10 28 cm 3 . 
If the annihilated field is limited to a smaller volume, the field 
strength must be even larger. The flare onset is quite rapid and 
the conversion of magnetic energy into thermal energy occurs 
within 20 - 40 s at most. 



3.2. LP 41 2-31 in the context of stellar flares 

The so far most detailed example of a coordinated optical and 
X-ray observation of a large flare on an M dwarf presented in 
the literatu re is an XMM-New ton observation of Proxima Cen 
(M5.5e) bv lGiidel etail fc004l) . The peak X-ray luminosities of 
the flares on Proxima Cen and LP 4 12-31 are quite similar and 
it is instructive to compare them in more detail. As in LP 412- 
31, the optical flare on Proxima Cen was shorter than the X-ray 
event, in line with the expectation from a flare scenario where 
the optical burst is considered as a proxy of the initial (impul- 
sive) heating of chromospheric plasma by accelerated electrons. 
However, there are marked phenomenological differences be- 
tween the two events: The Proxima Cen X-ray flare was a long- 
duration event with a rise time of 12 min and decay time of a sig- 
nificant fraction of a day, in contrast to the short flare observed 
on LP 4 12-31. The onset of the optical event preceeded that of 
the X-ray flare in the case of Proxima Cen while in LP 412- 
31 any time lags are below the time resolution of 20 s. In the 
chromospheric evaporation scenario stellar flares are expected 
to peak in the optical before the X-rays, but due to the lack of 
simultaneous multi band data observational evidence was foun d 
only in few cases dSchmitt et all 1 9931 iMitra-Kraev et a 1. 2005). 

The EPIC light curve of Proxima Cen displays - unlike the 
case of LP 412-31 - next to a major flare with total energy 
~ 1.5 ■ 10 32 erg, several small events with energies of a few 10 29 
to w 10 31 erg. Such observations of lower level flaring have lead 



to the idea of stellar coronal nano-flare heating. Studies of solar 
flares and some bright flare stars suggest the energy distribution 
of flares to foll ow a power law with slope « -2 in its differ en- 
tial form (e.g., Aschwa nden et alJl2000l lAudard et alJl2000t) . A 
systematic analysis of X-ray flaring in a sample of pre-main se- 
quence stars in the Taurus star forming region has yielded similar 
results ( Stelzer et al] 120061) . In the context of LP 412-31 arises 
the question if the observed giant flare can be interpreted as the 
"high energy t ail" of a canonic al flare frequency distribution. 
According to lStelzer et alJ J2006) we expect to be sensitive to X- 
ray events with relative amplitude A K \ = {R mdx - Rq)/Rq > 0.3, 
where R mdx and Rq are the maximum and the quiescent soft X- 
ray count rate, given the sensitivity limit imposed by Poisson 
statistics and the observed quiescent EPIC/pn count rate of 
LP 412-31. This limit corresponds to flares with peak luminos- 
ity of log Lx [erg s -1 ] ~ 27.3. Assuming a duration of 5000 s for 
such micro-flares - probably a conservative estimate consider- 
ing the short duration of the observed giant flare - the emitted 
energy of such events would be 3 ■ 10 30 erg. If the number dis- 
tribution of flares on LP 4 12-31 were as observed for other stars, 
one would expect to see » 100 such events for each giant flare, 
in clear contrast to our EPIC/pn data. Depending on the absolute 
flare frequency, it would in principle be possible that we failed 
to detect smaller flares because of the limited observing time, 
however, it is as well possible that the actual flare energy dis- 
tribution of LP 4 12-31 is rather flat. In this case, the origin of 
the quiescent X-ray emission (often considered to be composed 
of unresolved micro- and nano-flares) remains unclear. Future 
multi-wavelength monitoring observations of ultracool dwarfs 
will prove useful to constrain the structure and energetics of their 
outer atmospheres. 
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